METHOD AND ARRANGEMENT FOR MEDICAL X-RAY IMAGING 
BACKGROUND OF THE INVENTION 
[0001] Three-dimensional X-ray imaging is based on taking 

several 1-D or 2-D projection images of a 3-D body from different directions. If 1- 
5 D projection images are available from all around a 2-D slice of the body with 
dense angular sampling, the inner structure of the slice can be determined. This is 
known as Computerized Tomography (CT) imaging technology, which is widely 
used in medicine today. A crucial part of CT technology is the reconstruction 
algorithm taking the X-ray images as argument and producing a voxel 
10 representation of the 3-D body. 

[0002] In many practical cases X-ray projection images are 

available only from a limited angle of view. A collection of X-ray images of a 3-D 
body is called sparse projection data if (a) the images are taken from a limited angle 
of view or (b) there are only a small number of images. Sparse projection data does 
15 not contain sufficient information to completely describe the 3-D body. 

[0003] However, some a priori information about the body is 

typically available without X-ray imaging. Combining this information with sparse 
projection data enables more reliable 3-D reconstruction than is possible by using 
only the projection data. 
20 [0004] Traditional reconstruction algorithms such as filtered 

backprojection (FBP), Fourier reconstruction (FR) or algebraic reconstruction 
technique (ART) do not give satisfactory reconstructions from sparse projection 
data. Reasons for this include requirement for dense full-angle sampling of data, 
difficulty to use a priori information, for example nonnegativity of the X-ray 
25 attenuation coefficient, and poor robustness against measurement noise. For 
example the FBP method relies on summing up noise elements with fine sampling, 
leading to unnecessarily high radiation dose. 

[0005] The traditional setting of the image reconstruction 

problem in tomography assumes that the imaging geometry is known accurately. 
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The prior art requires either reference feature points, with which the unknown part 
of the imaging geometry (e.g. motion of the imaging device or the object) is 
computed, or precise prior knowledge of the imaging geometry. Anyway the prior 
art does not know effective methods to determine the imaging geometry like the 
5 motion of the imaging device and/or the object, though all parts of the imaging 
geometry are rarely precisely known in advance. 

BRIEF DESCRIPTION OF THE INVENTION 
[0006] The aim of the invention is to overcome the problems met 

10 in 3-D reconstruction of the object that occur when using traditional reconstruction 
algorithms with sparse projection data. This is achieved by a method for producing 
three-dimensional information of an object in medical X-ray imaging. The object 
is X-radiated from at least two different directions and the said X-radiation is 
detected to form projection data of the object. The object is modelled 

15 mathematically utilizing the projection data to solve the |[sl]imaging geometry|[s2] 
and/or the motion of the object. Said projection data and said mathematical 
modelling of the object are utilized in Bayesian inversion based on B ayes' formula 
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p{x,0 | m) = 



Ppr(0)P P r( x )P( m \ x >&) 

p(m) 



to produce three-dimensional information of the object, the prior distribution p pr (0) 
representing the prior knowledge of the imaging geometry] and/or the motion of the 
object, the prior distribution p pr (x) representing mathematical modelling of the 
object, x representing the object image vector, which comprises values of the X-ray 
25 attenuation coefficient inside the object, G representing the parameter vector of the 
imaging geometry| and/or the motion of the object, m representing projection data, 
the likelihood distribution p(m|x,0 ) representing the X-radiation attenuation model 
between the object image vector x, geometry parameter vector 9 and projection data 
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m, p(m) being a normalization constant and the posteriori distribution p(x,0|m) 
representing the three-dimensional information of the object and the imaging 
geometry including the motion of the object). 

[0007] The invention also relates to a medical X-ray device 

5 arrangement for producing three-dimensional information of an object (4) in a 
medical X-ray imaging, The medical X-ray device arrangement comprises: 

- an X-ray source for X-radiating the object from at least two different 

directions 

- a detector for detecting the X-radiation to form projection data of the 

10 object 

- means for modelling the object mathematically |[s6]utilizing the 
projection data to solve the imaging geometry and/or the motion of the object, 
where the solving concerns either some or all parts of the imaging geometry and/or 
the motion of the object. 

15 and the medical X-ray device arrangement comprises: 

- means for utilizing said projection data and said mathematical 
modelling of the object in Bayesian inversion based on Bayes' formula 
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p{x,6 | m) = 



^P P r(0)PprWP( m \ X >0) 



p(m) 



to produce three-dimensional information of the object, the prior distribution ppr(9) 
representing the prior knowledge of the the imaging geometry and/or the motion of 
the object (4), the prior distribution ppr(x) representing mathematical modelling of 
the object, x representing the object image vector, which comprises values of the 
25 X-ray attenuation coefficient inside the object, 0 representing the parameter vector 
of the imaging geometry and/or the motion of the object (4), m representing 
projection data, the likelihood distribution p(m|x,9 ) representing the X-radiation 
attenuation model between the object image vector x. geometry parameter vector 0 
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and projection data m, p(m) being a normalization constant and the posteriori 
distribution p(x,0|m) representing the three-dimensional information of the object 
(4) and the imaging geometry including the motion of the object. 

[0008] In this invention, the imaging geometry and/or the motion 

of the object are solved with the image reconstruction simultaneously or having 
different timings. This is based on that projection data itself simultaneously 
contains sufficient amount of information about the imaging geometry and the 3D 
reconstruction. Neither reference feature points nor prior knowledge of the motion 
are needed. 

[0009] The invention is also based on that biological issues have 

that kind of statistical a priori information that this information can be successfully 
utilized with Bayesian inversion in medical x-ray imaging. The suitable a priori 
information makes possible to model the biological tissue mathematically 
accurately enough and independently of X-ray imaging. From biological tissue it is 
possible to compile qualitative structural information which makes it possible to 
utilize the Bayesian method successfully to solve the problems in medical three- 
dimensional x-ray imaging. There is certain regularity in biological tissues and this 
regularity is useful especially with the Bayesian method. 

[0010] For example 10 x-ray images are taken from breasts of 

different persons. From these x-ray images it is noticed that there is much 
similarity in the statistical structure of the breasts between different people. In 
other words biological tissues and x-ray images taken from the biological tissues 
have similar or almost similar statistical structure between different persons. 

[0011] With Bayesian inversion it is possible to utilize a priori 

information efficiently in 3-D reconstruction from sparse projection data, because 
the suitable a priori information from the biological tissue makes possible to model 
the biological tissue mathematically accurately enough and independently of X-ray 
imaging. In addition, a priori information of the imaging geometry, such as 
approximate positions of the radiation source and the detector, can be similarly 
taken into account in a sound, statistical way. Any collection of projection data can 



4 



be used in the reconstruction. Application-dependent a priori knowledge can be 
used to regularize the ill-posed reconstruction problem. 

[0012] This invention improves the quality of 3-D reconstructions 

over traditional methods. In addition, the number of radiographs can be minimized 
without compromising quality of the reconstruction. This is very important in 
medical applications since the X-ray dose of the patient can be lowered. 



BRIEF DESCRIPTION OF THE FIGURES | 
[0013] Figure 1 shows a simple example of X-ray imaging. 

10 [0014] Figures 2 A - 2B show a difference between global 

tomography and local tomography. 

[0015] Figures 3 A - 3B show different types of sparse projection 

data. Every black dot represents a location of the X-ray source for taking one 
projection image. 

15 [0016] Figures 4 A - 4B show examples of parts of boundaries and 

cracks of an object visible and invisible without a priori information. 

[0017] Figure 5 illustrates "pencil beam" X-ray attenuation 

model. 

[0018] Figure 6 is basic flow chart of the method according to the 

20 invention. 

[0019] Figure 7 is an intraoral X-ray device arrangement 

presenting one preferred embodiment of the invention. 

[0020] Figure 8 shows a measurement geometry of dental 

limited-angle tomography with a digital intraoral sensor. 
25 [0021] Figure 9 shows results of 3D reconstruction of head 

phantom in the first preferred embodiment. 

DETAILED DESCRIPTION OF THE INVENTION 
[0022] In conventional tomography, the interior of an object is 

30 reconstructed from tomographic projections such as X-ray images. All the current 
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image reconstruction methods assume that projection, geometry of the imaging 
device is either known or solved in advance by using fiducial or non-fiducial 
feature points in the images. With this invention a novel approach is presented 
where the imaging geometry is solved with the reconstruction problem while no 
5 correspondence information is needed. This approach is a direct application of 
Bayesian inversion theory and produces the maximum likelihood or maximum a 
posteriori estimates for the imaging geometry like motion parameters under the 
selected noise and prior distributions. This invention can be utilized assuming 
either ID or 2D projections. To utilize the invention with ID projections is also 

10 often a sufficient assumption since the reconstruction problem is frequently 
performed as a stack of 2D reconstruction problems for computational reasons. 
However, the method can be also directly extended to the general 3D case. 

[0023] In practical imaging situations X-ray images are not 

always available from all around the body. The body might be visible only from 

15 certain directions due to imaging geometry. For example this is the case in 3-D 
mammography with the breast compressed against the detector, or in intraoral 
dental imaging with the detector inside the patient's mouth. This situation is called 
limited-angle tomography. Also, the region of interest might be surrounded by 
tissue that need not be imaged, like in extraoral dental imaging. This situation is 

20 called local tomography. In addition the number of radiographs should be 
minimized in medical applications for reducing the X-ray dose of the patient. 

[0024] In the preferred embodiments of the invention Bayesian 

inversion algorithms are used to create a new type of 3-D medical X-ray imaging. 
It is intermediate between a projection radiograph and a full computed tomography 

25 scan. Two steps are needed to perform successfully such imaging: In step one, the 
doctor (a) collects all the prior information he or she has on the tissue of interest, 
and (b) takes the minimum number of radiographs containing the desired 3D 
information. In step two, a tomographic algorithm based on Bayesian inversion is 
used for 3D reconstruction. The algorithm contains two separate mathematical 

30 models. First, all a priori information (i.e. information that is independent of X-ray 
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imaging) is used to model the unknown object mathematically. The model is put 
into the form of prior distribution in Bayes' formula. In addition, prior information 
of the imaging geometry parameters can be modelled. Second, the measurement is 
modelled mathematically. This involves the geometry of the imaging device, 
5 position of the detector and X-ray source during the exposure of each projection 
image and a model for X-ray attenuation inside tissue. The mathematical model of 
the measurement is put into the form of likelihood distribution in Bayes 5 formula. 

[0025] In figure 1 is shown a simple example of X-ray imaging, 

where an X-ray source 2 is placed on one side of an object 4 under imaging. 

10 Radiation passes through the object and is detected by a detector 6 on the other 
side. The X-ray source is for example a X-ray source of an intraoral X-ray source 
of a dentist, of a dental panoramic X-ray device, of a surgical C-arm X-ray device, 
of a mammography device or of any other medical X-ray device and the detector 6 
is a detector of some of those devices. Usually the detector 6 is a digital sensor that 

15 can be thought of as a 2-D array almost pointlike detectors. 

[0026] The 3-D body under imaging is modelled by nonnegative 

X-ray attenuation coefficient. The value gives the relative intensity loss of the X- 
ray travelling within a small distance dr: 



20 y=-Jt(r)<fr (1) 



[0027] The X-radiation has initial intensity I 0 when entering the 

object 4 and a smaller intensity I 0 when exiting the object. The following equation 
shows the attenuation law: 



25 



I 7^ dr=logI °~ l1 (2) 
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where initial intensity I 0 is known by calibration, and intensity after object I x is 
known from the corresponding point value in a projection image. Thus the 
measured data is the integral of x along the line L. 

[0028] In the above model it is not taken into account (a) 

5 scattering phenomena resulting in X-rays changing their direction, or (b) the 
dependency of attenuation on the X-ray spectrum resulting in low-energy photons 
being more easily attenuated than high-energy ones. Effect (b) causes errors in 
measurements and is sometimes referred to as beam hardening. More detailed 
models can be utilized in Bayesian inversion to overcome (a) and (b). 

10 [0029] In medical imaging the geometrical arrangements of the 

X-ray source and digital sensor vary according to the diagnostic task and 
equipment. Figures 2A-2B and 3 illustrate the types of tomographic data resulting 
from different imaging situations. For clarity, here are presented two-dimensional 
examples; similar situations can be considered in 3-D. 

15 [0030] In figure 2A-2B are shown two cases according to whether 

the whole object 4 is fully visible in each projection or not. These cases are called 
global tomography and local tomography, respectively. In figure 2B is radiated 
only ROI (Region Of Interest) 8. 

[0031] The types of data described in figures 3 A - 3B cover a 

20 large range of specific imaging tasks. The choice of data collection dictates what 
kind of features and details it is easiest to reconstruct reliably from the data. Thus 
the choice of projection views must be made according to the diagnostic task at 
hand. In traditional CT imaging, projections are taken from all around the object 4. 
In the preferred embodiments of the invention radiation dose is lowered by 

25 sampling the angular variable more sparsely. Figure 3A presents more dense 
angular data with full angle and limited angle measurements. Figure 3B presents 
more sparse angular data, also with full and limited angles of measurement. 

[0032] In figures 4A and 4B there are presented examples of parts 

of boundary and cracks that are visible or indetectable in reconstruction without a 

30 priori information. Figure 4A there is an object 4 under imaging with edge on the 
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surface of the object. The edge is detectable on the leftmost position where the 
edge is more parallel to the direction of the X-rays. In the figure 4B there is a crack 
inside of the object 4. The crack is detectable in the leftmost position parallel to the 
angle of measurement. 
5 [0033] The main idea in Bayesian inversion method is to consider 

the inverse problem as a problem of statistical inference. All variables are 
redefined to be random variables. The randomness reflects uncertainty of their 
actual values and the degree of uncertainty is coded in the probability distributions 
of these random variables. 
10 [0034] Typically, the model for X-ray tomography problem is 

represented by the formula 

m=Ax+e (3) 

15 where m is a vector of projection measurements (all the projection images 
represented as one vector), x is a vector that represents the object densities (voxel 
values) at given positions in 3D space, A is a transformation matrix that 
characterises how the projections are constructed from the 3D object, hence 
contains the imaging geometry, and e is a random, noise vector. The variables m, x 

20 and e are random variables and the matrix A comes from the pencil beam model for 
the X-ray attenuation: A contains the weights describing how much each voxel 
contributes to the X-ray attenuation along the X-ray paths. In this way the integral 
in formula (2) is approximated. This is illustrated by an example in figure 5. 

[0035] In the traditional setting, the objective is to reconstruct the 

25 structure (density) of the object given the projection images and the imaging 
geometry. As an extension to the traditional setting, in this invention both the 
reconstruction and the imaging geometry are solved simultaneously. To emphasise 
that, we write the model (3) in the form 

30 m = A e x + e, 
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where we just explicitly say that A depends on the geometry parameters G that 
define the unknown part of the imaging geometry and motion of the object. Hence, 
the unknowns in the equation are both 0 and x. 
5 [0036] Simultaneous solving of both 0 and x is possible, because, 

the model (3) constrains the imaging geometry sufficiently. This is due that, even 
though there are more unknowns than equations, the equations are not independent 
and, if e.g. the motion parameters are incorrect, the equations are also inconsistent. 
Moreover, if the geometry parameters 0 were incorrectly fixed, in general, there is 

10 no reconstruction x that would produce the observed images m. Therefore, 
excluding certain special cases, the geometry parameters are completely determined 
from the noise model, though the image object x is not uniquely (additional prior 
information is required), as far as the imaged object is sufficiently complex. 
Therefore any feasible reconstruction x (in the likelihood sense) can be used to 

15 compute 0, after which the reconstruction x can be updated to accord best with the 
prior information. 

[0037] Assume now that the object image vector x, geometry 

parameter vector 0, and the noise are jointly independent random variables. The 
conditional probability distribution of x given the measurement m is given by 
20 Bayes' formula in the form 



p{x,0 | m) = 



_P^O)P^{x)p{m\x,0) 



p(m) 



(4) 



where p(m) is a normalization constant. The density p(x,0|m) is called the 
25 posteriori distribution of x and 0. p(m|x,0) is a likelihood distribution representing 
the X-radiation attenuation model between the object image vector x, geometry 
parameter vector 0, and projection data m. 
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[0038] The density P pr (x), called the prior distribution of x, is 

designed to contain all possible information available of the actual object 4 
independently of X-ray imaging. The density P pr (0) is the prior distribution of 0, 
and it models the available information of the imaging geometry including the 
5 motion of the object. If we do not have any such prior information of the 
parameters 0, we may use the non-informative prior P pr (6) = constant. It is crucial, 
in contrast to several classical regularization methods, that the choice of the prior 
distributions should not be based on the projection data. The rule of thumb in the 
design of prior distributions is that typical image vectors (say, of some existing 

10 library) should have high prior probability while atypical or impossible ones should 
have low or negligible probability. 

[0039] In the framework of the Bayesian inversion theory, the 

posterior distribution in formula (4) represents the complete solution of the 3D 
reconstruction problem as well as of the imaging geometry estimation problem. 

15 [0040] To produce an image of the object 4, and similarly an 

estimate for the geometry parameters, based on the posterior distribution, several 
alternatives exist. The most common ones are the maximum a posteriori estimator 
(MAP) and conditional mean estimator (CM). Let y = (x,0) i.e. y is the vector 
composed of both the object image and geometry parameter vectors. The 

20 estimators are defined by the formulas 

XVmap I m ) = max/?(> | m) (5) 

where the maximum on the right hand side is taken over all y, and 

25 

^cm = \yP (y I m )dy 

[0041] Finding the MAP estimator is an optimization problem 

while finding the CM estimator is a problem of integration. 
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[0042] In figure 6 is presented a basic flow chart of the method 

according to the invention. In method step 600 the object is modelled 
mathematically independently of X-ray imaging. In method step 602 the object is 
5 X-radiated from at least two different directions. In method step 604 the said X- 
radiation is detected to form projection data of the object. In method step 606 said 
projection data and said mathematical modelling of the object are utilized in 
Bayesian inversion based on Bayes formula to produce three-dimensional 
information of the object. The method step 600 is also possible to perform after 
10 method step 602 or 604. 

[0043] In the first preferred embodiment of the invention is 

presented an application to dental radiology. 

[0044] X-ray projection images are conventionally used in dental 

radiology. However, certain diagnostic tasks require more precise knowledge of 
15 the 3D structure of tissue than is available in two-dimensional radiographs. Such 
tasks include implant planning and detection of bone loss between tooth roots. 

[0045] In figure 7 is presented an intraoral X-ray device 5 

arrangement presenting the first preferred embodiment of the invention. It is 
important to notice that this is only an example of the medical X-ray device 5 
20 arrangement where the invention is possible to be utilized. 

[0046] The medical x-ray device 5 in the preferred embodiments 

of the invention is for example a dental panoramic X-ray device, a surgical C-arm 
X-ray device or a mammography device. 

[0047] In the first preferred embodiment of the invention the 

25 articulated arm arrangement 3 moves the X-ray source 2 to the right position. The 
X-radiation begins by pressing the exposure button 12. The X-ray source 2 X- 
radiates the object 4, which is for example teeth of a patient. The detector 6 detects 
the X-radiation. The image information which is got by detecting the X-radiation is 
sent by communication link 16 to the computer 14. The computer comprises the 
30 software means 15 to process the image information according to the invention. 
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There can be more than one computer 14 and also the software means 15 can 
situate in more than one computer 14. For example the first computer 14 is 
computer which is used in x-ray imaging. The second computer 14 is computer 
which is used in processing the image information according to the invention. It is 
5 possible to have the second computer 14 far away from the actual medical x-ray 
device 5. For simplicity in figure 7 is shown only one computer 14. 

[0048] In the first preferred embodiment (Fig. 7) of the invention, 

the dentist's X-ray equipment is used for taking a set of 2D projection images that 
are used as input for Bayesian 3D reconstruction algorithm. Such equipment 
10 includes an intraoral X-ray unit and a digital intraoral sensor. 

[0049] Benefits of this approach over conventional CT scan are 

low cost and convenient usage, 
- high resolution of the projection images, 

possibility to take as few radiographs as is needed to capture the 
15 relevant 3D information, minimizing radiation dose to the patient, 

the possibility to choose imaging directions so that the X-rays do 
not pass through the whole head but only the interesting tissue, further reducing 
dose. 

[0050] In the preferred embodiments of the invention, Bayesian 

20 inversion is used for the 3D reconstruction. Input data for the algorithm is the set of 
projection images and the following a priori knowledge: 

(a) Dental tissue consists of few approximately homogeneous 
regions with well-defined boundaries. 

(b) Dental tissue can only attenuate X-radiation, not intensify 

25 it. 

[0051] The example of the detector 6 used in the first preferred 

embodiment of the invention is based on charge coupled device (CCD) technology 
and has dynamic range of 4096 gray levels. The size of the active imaging area is 
34 mm * 26 mm and the resolution is 872 * 664 pixels. After exposure, each pixel 
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contains an integer which is proportional to the number of X-ray quanta that hit the 
pixel's area. 

[0052] Alternative detectors include any other digital intraoral 

sensor, digitized X-ray film, or any intraoral sensing device converting detected X- 
5 ray photons to a digital image. 

[0053] In the measurement geometry the focal point of the X-ray 

source 2 moves horizontally on a circle in the plane 10, with center at the detector 
6, see Figure 8. In Figure 8, the detector 6 corresponds to one row of CCD pixels 
in the CCD sensor. 

10 [0054] In figure 8 is presented an example of taking seven 

radiographs of some teeth of a patient. The total opening angle of the projection 
views is 55 degrees in this example. Figure 8 represents the xy plane, and the z- 
coordinate axis is thought to be perpendicular to the plane of the paper. 

[0055] In the preferred embodiments, mathematical modelling of 

15 the object 4 can be done with a prior distribution of the form in formula (7): 

P p Xx)=^V^YPn(x)) (7) 

N 

where the sum is taken over a collection of 3D neighbourhoods N and the value 
20 U N (x) depends only on the values of voxels belonging to the neighborhood N, and 
a is a positive regularization parameter used to tune the width of the prior 
distribution. The neighborhood N is typically a set of voxels whose center points 
are closer to each other than some predefined maximum distance. 

[0056] It is to be noted that formula (7) does not define a 

25 probability distribution since it does not integrate to 1. Instead, the integral 
diverges. However, when used in the Bayes formula (4) with a suitable likelihood 
distribution, the posterior distribution is integrable. The same remark concerns the 
priors derived from (7), that is, formulae (9), (10) and (11). 
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[0057] In the preferred embodiments of the invention the three- 

dimensional problem can be reduced to a stack of two-dimensional problems each 
corresponding to a plane determined by a constant value of z. Here figure 8 
represents exactly the situation in the xy plane, i.e. the plane z=0. Each row in the 
5 detector corresponds to one such 2D problem. This approach leads to small 
approximation error because the X-ray source does not move in the correct plane 
for nonzero z, and this error is neglected. 

[0058] Next the modeling of a 2D tomographic problem is 

explained. Let 

10 

m(j) = A xG) + eG) (8) 

denote the jth 2D tomographic problem. Here the vector m(j) contains the readings 
on jth row from each of the seven radiographs. The vector x(j) is the jth slice of the 

15 3-D representation x of the object 4 under imaging. The matrix A is the pencil 
beam model matrix for the 2-D tomography problem. Now x(j) is a 2D array of 
pixels. The pixels are denoted by x(j)[k,q], where k=l,2,...,K is row index and 
q=l,2,...,Q is column index. 

[0059] The mathematical modelling of the object 4, i.e. 

20 incorporation of prior information, is next explained for the 2D slice. Define 

P pr (^ 0) )=exp^E fu(N(xV% q ]))) 

k=l q=\ 

where the value U(N(x(j)[k,q])) depends only on the values of pixels belonging to 
25 the neighborhood N of the pixel x(j)[k,q] and a is a positive regularization 
parameter used to tune the width of the prior distribution. The choice of the 
neighborhoods N(x(j)[k,q]) is arbitrary, but a typical choice is the set of pixels 
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whose centerpoints are closer to the centerpoint of x(j)[k,q] than some predefined 
distance. 

[0060] In the preferred embodiment, the neighborhoods are 

chosen to consist of two adjacent pixels. Further, the function U is chosen to be a 
5 power of absolute value of the difference between adjacent pixels. These choices 
lead to formula (10): 



P pr (x^)=exp 



K-\ Q K Q 



-<£ X \x U %q]-x"[k + l, g y + + X X \x U) [k,q]-x"[k, g+ l]*) 



k=\ q=\ k=\ q=] 



(10) 



10 where s is a positive real number and a is a regularization parameter used to tune 

the width of the prior distribution. 

[0061] In the equation (11) is presented non-normalized total 

variation (TV) distribution, when s=l in equation (10). This is the case according 

to the preferred embodiment of the invention, because when s=l the prior model 
15 becomes such that it gives high probability for objects consisting of a few areas of 

different attenuation with well-defined boundaries. 



prTJix u) )=ex$ 



f -^L t \x°%g]-^[k^g] ++± Z \x u % q ]-^[k 9 q+l] )' 

^ *=1 q=\ Ar=l q=l ) 



(11) 



20 X-radiation can only attenuate, not strengthen, inside tissue. This leads to 
positivity prior pos defined by 

pos(x(j)) = 1 if all pixels of x(j) are positive, 0 otherwise (12) 

25 [0062] The explanation of modeling of one 2D tomographic 

problem is now complete. 

[0063] The stack of 2D problems are connected to each other by 

demanding that consecutive 2D slices x(j) and x(j-l) should not be very different. 
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[0064] Mathematically this is expressed as follows: 

pr3D(xG)) = exp(- y | xG)[k,q] - xG-l)[k,q] |), (13) 

where the sums are taken over k=l,...,K and q=l,...,Q and y > 0 is another 
regularization parameter. So the prior distribution in formula (4) is the product of 
(11), (12) and (13): 

Ppr(xG)) = prTV(xG))pos(xG)) P r3D(xG)). (13) 

[0065] The measurements are taken into account in the form of 

the likelihood distribution 

P(m 0) |w°Vcexp{-^ (m U) -Ax U) )} (15) 

where S is the covariance matrix of the Gaussian noise vector e. c ia a 
normalization constant. 

[0066] In reality, X-radiation measurement noise is Poisson 

distributed. Formula 15 uses a Gaussian approximation to the Poisson distribution. 

[0067] All the parts are presented of the right hand side of (4). 

Next it is looked for the MAP estimate of xQ), that is, the image xG) that gives the 
largest value for the posteriori distribution (4). This is equivalent to finding the 
image xG) that minimizes the expression 

F(x U) =^(m U) -Ax U) )^ ( miJ) -AX U) + -log(prTV(x U) ))-log(pr3D(x U) )) (16) 
with the additional requirement that every pixel of xG) is positive. 
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[0068] Minimization of (16) is difficult since F is not a 

differentiate function due to the absolute values present in (11) and (13) and due 
to the sharp cutting in (12). The absolute values in (1 1) and (13) are replaced by 

5 M * hp (0 = log(cosh(^)) ( j 7) 

where (3>0, enabling the use of efficient gradient-based minimization routines. The 
positivity constraint is taken care of by an interior search minimization method. 

[0069] The results are shown in figure 9 in comparison with 

10 traditional tomosynthesis technique based on backprojection. 

[0070] In this description explained mathematic modelling and 

other image information processing are performed by software means 15 to 
accomplish three-dimensional information of the object 4. 

[0071] In the second preferred embodiment the invention is 

15 utilized in dental panoramic x-ray imaging. The object 4 is typically teeth of the 
patient and the medical x-ray device 5 arrangement is dental panoramic x-ray 
device arrangement. 

[0072] In the third preferred embodiment the invention is utilized 

in mammography, for example in FFDM (Full Field Digital Mammography). 
20 There the object 4 is breast of human being and the medical x-ray device 5 
arrangement is mammography device arrangement. 

[0073] In the fourth preferred embodiment the invention is 

utilized in surgical x-ray imaging. There the object 4 is the body part under surgery 
and the medical x-ray device 5 arrangement is surgical x-ray device arrangement. 
25 [0074] In the first, second, third and fourth preferred 

embodiments of the invention the basic method steps are same as mentioned with 
the flow chart in figure 6. The utilizing of the invention in the second, third and 
fourth preferred embodiment is similar to what is described with the first preferred 
embodiment of the invention and elsewhere in this application except different 
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medical x-ray imaging applications and their differences because of different 
medical x-ray devices and different objects to be x-ray imaged. 

[0075] What makes this invention especially inventive is that the 

imaging geometry and/or the motion of the object is solved with the image 
5 reconstruction simultaneously or having different timings. This is based on that 
projection data itself contains sufficiently information about the imaging geometry 
and also sufficiently information about the 3D reconstruction. Neither reference 
feature points nor prior knowledge of the motion are needed. 

[0076] This can be done as described above with the preferred 

10 embodiments and with the means 15 for modelling the object 4 mathematically 
utilizing the projection data to solve the imaging geometry and/or the motion of the 
object and with the means 15 for utilizing said projection data and said 
mathematical modelling of the object in Bayesian inversion based on Bayes' 
formula 
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p(x,0 | m) = 



_P P r(0)P P r( X )P( m \ X >0) 



p(m) 



to produce three-dimensional information of the object, the prior distribution p pr (0) 
representing the prior knowledge of the imaging geometry and/or the motion of the 

20 object 4, the prior distribution p pr (x) representing mathematical modelling of the 
object, x representing the object image vector, which comprises values of the X-ray 
attenuation coefficient inside the object, 9 representing the parameter vector of the 
imaging geometry and/or the motion of the object, m representing projection data, 
the likelihood distribution p(m|x,9 ) representing the X-radiation attenuation model 

25 between the object image vector x and projection data m, geometry parameter 
vector 0, p(m) being a normalization constant and the posteriori distribution 
p(x,0|m) representing the three-dimensional information of the object 4 and the 
imaging geometry including the motion of the object. 
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[0077] Although the invention is described above with reference 

to the examples illustrated -in the attached figures, it is to be understood that the 
invention is not limited thereto but can be varied in many ways within the inventive 
idea disclosed in the attached claims. 
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